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GLOSSARY  OF  SYMBOLS  AND  ABBREVIATIONS 
a  Radius  of  the  earth  (-  6.37  megameters) 

c  Velocity  of  light  in  free  space  (•  3  x  108  meters  per  second) 

c/v  Earth-ionosphere  waveguide  phase-velocity  ratio 

ELF  Extremely  low  frequency  (30  to  300  Hz) 

Ep  Horizontal  electric-field  component  in  the  p  direction  (volts/meter) 

E^  Horizontal  electric-field  component  in  the  <j>  direction  (volts/meter) 

Ez  Vertical  electric-field  component  (volts/meter) 


f(x) 

1  -  xH0(2)(x)/H1(2)(x) 

G(u) 

(x)coth  u  *  i1  -  !)u2 

csch2  u 

G(t) 

(f)c°th  t .  (i  - 

csch2  t 

h 

Hq(2)(x) 

Ionospheric  reflection  height  (meters) 

Hankel  function  of  the  second  kind,  order  zero,  and  argument  x 

h!-:  1  (x) 

Hankel  function  of  the 

second  kind, 

order  one,  and  argument  x 

HED  Horizontal  electric  dipole 

HMD  Horizontal  magnetic  dipole 

H„  Horizontal  magnetic-field  component  in  the 

o  direction  (amperes/meter) 

Horizontal  magnetic-field  component  in  the 
?  direction  (amperes/meter) 

li:  t  i  Git)  *  V  ( t  ) 

I  Current  (amperes) 

k  2t/\q  (meters'1) 

m  Magnetic  dipole  moment  Tampere  meters-) 


NIKS 


Met er-k i logram- second 
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NUSC 

P 

S 

ikS0 

t 

T 

u 

V 

VED 

W 

x 

z 

'e 

Yi 

£o 

Ee 


Naval  Underwater  Systems  Center 
Electric  dipole  moment  (ampere  meters) 

,  spherical  earth  spreading  factor 

ELF  earth- ionosphere  waveguide  propagation  constant 


u  _  ttp 
(c/v)2  2h(c/v)2 


TTP 

2h 

Earth- ionosphere  waveguide  propagation  velocity  (meters/second) 
Vertical  electric  dipole 

££[»<«  *^C(u)x2H0(2)(x)] 

kp (c/v) 

Vertical  distance  in  a  cylindrical  coordinate  system  (meters) 
Earth-ionosphere  waveguide  attenuation  rate  (nepers/meter) 
ik,  free-space  propagation  constant  (meters':) 

[iuu0(oe  +  iwCg)]1^2  ,  propagation  constant  in  the  earth  (meters" *) 

[iwu0(CT^  +  iue^)]1^2  ,  propagation  constant  in  the  iono¬ 
sphere  (meters'1) 

10'9/36tt  farads/meter,  permittivity  of  free  space 
Effective  permittivity  of  the  earth  (farads/meter) 

Effective  permittivity  of  the  ionosphere  (farads/meter) 

120",  impedance  of  free  space  (ohms) 

,  impedance  of  the  earth  (ohms) 


e 


Free-space  wavelength  (meters) 
Wavelength  in  the  earth  (meters) 
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=  4ir  x  10"7  henries/meter,  permeability  of  free  space 

Radial  distance  in  a  cylindrical  coordinate  system  (meters) 

Distance  where  the  minimum  value  of  the  VED  radial  wave  impedance 
occurs  (kilometers) 

Effective  conductivity  of  the  earth  (Siemens/meter) 

Effective  conductivity  of  the  ionosphere  (Siemens/meter) 

Azimuth  angle  in  a  cylindrical  coordinate  system 
2rrf  radians/ second,  angular  frequency 


A  COMPILATION  OF  FIELD-STRENGTH  FORMULAS  FOR  ELF  RADIO-WAVE 
PROPAGATION  IN  THE  EARTH- IONOSPHERE  WAVEGUIDE 

INTRODUCTION 


It  is  the  purpose  of  this  report  to  present  new  formulas  for  horizontal 
electric  dipole  (HED),  horizontal  magnetic  dipole  (HMD),  and  vertical  electric 
dipole  (VED)  extremely  low  frequency  (ELF)  radio-wave  propagation  in  the  earth 
ionosphere  waveguide.  These  new  formulas  extend  the  results  of  Wait1-  and 
Galejs,^  which  are  valid  for  measurement  distances,  p,  greater  than  approxi¬ 
mately  three  ionospheric  reflecting  heights,  h,  down  to  the  quasi-nearfield 
range,  which  is  defined  as  the  range  where  p  is  greater  than  an  earth  wave¬ 
length,  \e,  but  much  less  than  a  free-space  wavelength,  X0.  For  the  sake  of 
completeness,  the  abovementioned  previously  derived  formulas  also  will  be 
included. 


The  three  dipole  antennas  (VED,  HED,  and  HMD)  are  situated  at  zero  height 
with  respect  to  a  cylindrical  coordinate  system  (p,p,z)  and  are  assumed  to 
carry  a  constant  current,  I.  The  axes  of  the  VED  and  HED  (of  dipole  moment  p) 
are  oriented  in  the  z  and  x  directions,  respectively,  while  the  axis  of  the 
HMD  (of  dipole  moment  m)  is  oriented  in  the  y  direction.  The  ionosphere  is 
located  at  height  z  >.  h,  while  the  earth  is  located  at  height  z  <_  0.  The  prop 
agation  constant  in  the  air  is  denoted  by 

Y0  [■  ik  =  i2ir/\0 ]  , 

whereas  the  propagation  constants  in  the  earth  and  ionosphere  are  denoted  by 
■e  [=  V  3(:e  +  ] 

and 

Yi  [=  V iwu0  Oi  ♦  iwei)j  , 

respectively.  The  magnetic  permeability  of  the  whole  space  is  assumed  to 
equal  Uq,  the  permeability  of  free  space.  Meter-kilogram-second  (MKS)  units 
are  employed  and  a  suppressed  time  factor  of  exp(iujt)  is  assumed. 


DERIVATION  PROCEDURE 


Accounting  for  ionospheric  reflection  effects  out  to  distances  of  approx¬ 
imately  three  ionospheric  reflecting  heights,  h,  is  a  tedious  process  involv¬ 
ing  an  infinite  sum  of  images. However,  by  following  the  procedure 
outlined  by  Martin^  and  Bannister  and  Williams,11  we  find  that  each  VED,  HED, 
and  HMD  field-component  expression  can  be  multiplied  by  one  of  the  following 
four  functions: 


T 


V*  V* 


7^  TV  1^  7?  7 v  ■„"  ^  ^  >■  r-  ^  T^V"V^  T  w  '.' 


'.’ •  ''j  w"- v-- rwir: v 


T^7 


v  ’v  -jj  »  ;j  •  pp? 
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G(u) 

=  0^coth  u  +  ^1  -  “■ju2  csch2  u 

(1) 

G(t) 

=  0p-jcoth  t  +  ^1  -  jjt2  csch2  t 

(2) 

V(t) 

=  t2  coth  t  csch2  t 

(3) 

H  (t) 

*  G (t)  +  V(t)  , 

(4) 

where 


u 


TT£_ 

2h 


and 


t  =  - T  , 

(c/v)- 

where  c/v  is  the  earth-ionosphere  waveguide  phase-velocity  ratio 
c  -  3  x  108  m/s  . 


(5) 


(6) 


The  functions  G(t),  H(t),  and  V(t)  are  plotted  versus  t  in  figure  1.* 
Note  that  the  plot  of  G(t)  versus  t  is  also  a  plot  of  G(u)  versus  u.  When 
u  <  0.5,  G(u)  -  1;  when  t  <  0.5,  G(t)  ~  V(t)  -  1.  However,  when  t  <  2, 

Hft)  -  2.  Furthermore,  when  u  >  2.5,  G(u)  -  2u/-r;  when  t  >  2.5,  G(t)  -  2t/~. 
However,  when  t  >  4.5,  V(t)  ~  0  and  H(t)  -  2t/m. 


When  the  measurement  distance,  c ,  is  greater  than  approximately  three 

ionospheric  reflecting  heights  from  the  source,  each  VED,  HED,  and  HMD  field- 
component  expression  varies  as  a  Hankel  function, 


or 


(kSQo) 


or  3  combination  of  the  two),  where  k  =  2-/‘>-  and  iLS-  is  the  propagation 
constant  in  the  earth- ionosphere  waveguide.  is  related  to  the  phase 

velocity  v  and  attenuation  rate  a  by  the  formulas  c/v  =  RcSq  and  t  =  -S.~k  ImS 

At  ELF  (i.e.,  30  to  500  Hz),  Re(kS0c)  >>  Im(kS0o).  Therefore, 

if!,-  1  (kSr  1  -  11*  :j  [ \)e'  xz  (  '  ! 


All  figures  have  been  placed  together  at  the  end  of  this  report. 
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and 

H[2)(kS0p)  -  Hf2)(x)e‘ap 


(8) 


where 


x  =  kp(c/v)  .  (9) 

In  this  report,  we  will  use  previously  derived^-®  quasi-nearfield  range 
formulas  (p  >  Xe,  o  <<  XQ  ,  and  p  <  h/3),  as  well  as  the  Wait1  and  Galejs2 
formulas  (p  >  3h) ,  to  find  VED,  HED,  and  HMD  formulas  valid  at  ELF  for  p  >  Xg 
with  no  restrictions  on  the  ratio  of  p  to  h, 


As  an  example  of  our  derivation  procedure,  consider  the  VED  component. 
When  p  >  3h, 


hVE  ,  £g^)(4p)Hf2>  (X)  . 

For  x  <_  0.25,  equation  (10)  reduces  to 
„VE  .  PM  . 

^  2  ftp  2  \h  / 

When  p  <  h/3,  the  quasi -nearfield  range  formula  is 

hVE  .  P 

i>  2-tp2  * 

Since  G  ( u  J  -  1  for  c  <  h/3  and  c/h  for  c  >  3h ,  then,  for  x  <_  0.23. 

..VP  pGi'u) 


(10) 


(H) 


(12) 


Because  the  range  of  validity  of  equations  (10)  and  (13)  overlap  when 
c  >  3h  and  x  <_0.25,  re  can  substitute  G(u)  for  c/h  in  equation  (10)  to  obtain 
the  general  formula  valid  for  c  >  \e  (with  no  restrictions  on  the  ratio  of  c 
to  h) .  It  is 


pG  (u)  e"a0  /* i ’r xV ,  ( 2 )  f  -j 

(  1 

2-tc2  \  2  /H1  1  1  • 

^  1 

Wait's  expression!  for  H;  is 


Therefore , 

which ,  for  x  0 .25,  reduces  to 


'.VxV-'. 
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T  -  £(U I 
ikp 


(17) 


The  magnitude  of  T  (from  equation  (16))  is  compared  with  Wait's  infinite 
sum-of-images  result1  in  figure  2  as  a  function  of  distance  and  frequency. 

For  this  particular  comparison,  h  =  90  km  and  =  <®  (i.e.,  c/v  =  1.0  and 
a  =  0.0).  Note  that  the  agreement  is  excellent. 


We  know  from  previous  results1 >2  that  there  is  a  substantial  amplitude 
dip  in  the  VED  Ez  component  in  the  range  of  100  to  300  km  (depending  on  fre¬ 
quency)  .  Therefore,  we  will  let 


pVE  _  p  _  p 
Ez  czi  CZ2  * 


(18) 


For  p  <  h/3,  the  quasi-nearfield  range  formula  is 
„  ip  _  i60p 


2i\tae.Qpi  ko  ■ 

For  p  comparable  to  h,  we  can  easily  show  that 

F  i60pV (t) 

-1  '  ko3  ’ 


(!9) 


(20) 


which  reduces  to  equation  (19)  when  t  <  0.5  and  vanishes  when  t  >  4.5, 
When  o  >  3h, 

( c/v)  •-H  ( 3  ^  (x)  e"010  =  ^ i u x " H ;  ^ ^  (x)e~acj  . 

.25.  equation  21  i  reduces  to 


E. 


21 


ibOp 


(2u\  x2 ! 

in  f1-123^  ilT 

1 

i60pG (u) x2  | 

in  (lA2A  -  ±1 

\  x  /  2 

| 

*  kp 3  j 

V  x  >  2 

i 

m 

1 

(22) 


Therefore , 

E:.  -  . 

Employing  equations  (18),  (20),  and  (23)  results  in 


kc  :  I 


Vi  t 


1 1  u  )  x  H 


O'  lX]  ’ 


which  is  the  final  result. 

-  V  F  1 

Wait's  expression  for  EiL  is 
,  Vi:  IbOpk 


A  . 


J 
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Therefore, 

w  -  ^pr[v(t)  ♦  ^-G(u)x2h{2)(x)]  ,  (26) 

which,  for  x  <_  0.25,  reduces  to 

(CfH  -<*o*i[*.(Hr)-£]}.  (27) 

The  magnitude  of  W  (from  equation  (26))  is  compared  with  Wait's  infinite 
sum-of-images  result^  in  figure  3  as  a  function  of  distance  and  frequency. 

For  this  particular  comparison,  h  =  90  km  and  a^  =  00  (i.e.,  c/v  =  1.0  and 
a  =  0.0).  Note  that  the  agreement  is  excellent. 

From  equation  (27),  we  see  that  the  minimum  value  of  W  will  occur  when 
V(t)  -  G(u)x2£n(l .123/x) .  That  is, 

WMIN  '  X  G(uKc/v;)2  *  (28} 

When  u  >  2.5, 

WMIN  '  iu(c/v)2  =  i(^)(c/v)2  •  C29) 


FIELD-STRENGTH  FORMULAS 


In  this  section,  we  will  present  new  formulas  for  HED,  HMD,  and  VED  ELF 
radio-wave  propagation  in  the  earth-ionosphere  waveguide.  All  of  these  for¬ 
mulas  have  been  obtained  by  following  the  procedure  outlined  in  the  previous 
section.  They  are  valid  for  o  >  Ae,  with  no  restrictions  on  the  ratio  of  „ 
to  h. 


It  should  be  noted  that  for  2  Mm  <_  p  <_  19  Mm,  all  of  the  field-strength- 
component  formulas  presented  in  this  report  should  be  multiplied  by  the  spher¬ 
ical  earth  spreading  factor  S,  which  is  equal  to 


S 


o/a  1v 
sin(o/a) J  ’ 


(50) 


where  .  =  1/2  for  all  E, ,  Ez,  and  H.  components;  v  =  3/2  for  all  E.  and  H . 
components;  and  a  is  the  radius  of  the  earth  (~  6.37  Mm). 


FOR  THE  VED 


Expressions  for  the  VED  are 


nepG(u)e"-1D 

O  -  W 


l  a  1  ) 
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UHE 


L  PG (t)  cos  ^e-aoir/-i^x\u(2)f^gr^1 
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f(x) 


H0(2)  (x) 
h{2) (x) 


(40) 


The  magnitude  of  f(x)  is  plotted  in  figure  4  versus  x.  For  small  values 
of  x,  f(x)  -  1,  while  for  large  values  of  x,  f(x)  -  x. 


The  various  HED  components  are  related  by 


FHE 

uHE 

HP 

pHE 

P 

‘  7me  '  T)e  ’ 

(41) 

fHE 

c2 

in0xG(u) 

(42) 

uHE 

(c/v)f (x)G(t)  * 

pHE 

cz 

in0xG(u) 

(43) 

pHE 

cp 

ne(c/v)f (x)G(t)  ’ 

Hf 
hhh  ' 

E“E  -  -  f(x>  cot 

EHE  H(t)  ^  (p/a)  J 

(44) 

FOR  THE  HMD 

Expressions  for  the  HMD  are 


and 


(x)f(x) 


E»M 


mHM  .  mG  ( tj_  cos  9,0=  j  [ ’fzlZiy  [ ' 1  -  x  )  f  ( x )] 


(45) 

(46) 

(47) 

(48) 


(40) 
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The  various  HMD  field  components  are  related  by 

PHM  cHM 

H  EP 

hhm  '  "  hhm  '  ne  ’ 


E™1  in0xG(u) 

^HM  "  "  (c/v)f (x)G(t)  ’ 

4> 

E^1  in0xG(u) 

cHM  ne(c/v)f (x)G(t)  ’ 

hHM  pHM 

- -9H1  f(x)  cot  JLin(p/a)  1  . 

uHM  FHM  H(t)  rLXJ  c  z  9[  (p/a)  J 

no  H 

FIELD-STRENGTH  FORMULAS  FOR  x  <  0.25 


When  x  =  kp(c/v)  <_  0.25,  f(x)  -  1.0,  and  the  Hankel  functions  can  be 
approximated  by 


(^)h« 


(x)  -  1.0 


Hq  j (x) 


-  f  ]  • 


Thus,  the  VED  expressions  listed  in  the  previous  section  reduce  to 


eve  - 

nepG(u)e-ap 

to 
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pVE  oe(kp)G(u) 

CP 


(vCt)  -  gcu)4„(4“)  -  f])  . 


When  x  £  0.25,  the  HED  expressions  listed  in  the  previous  section  reduce 


eHE  ,  PG  (t)  cos  pe~ap 
p  2ir(ae  +  ituEgjp3 

eHE  _  _  pH(t)  sin  pe~ap 
^  2ir(a  +  iwee)p3 

ia)u0pG(u)  cos  pe'ap 

pnfc  _ 

2  "  2»Yep2 

hHE  _  pH(t)  sin  3e~ap 
p  2ttysP3 

hHE  _  pG(t)  cos  4>e~ap 

*  2iTYep 3 

E™  innkpG(u) 


E?c  ingkcG(u) 
pHE  '  neG(t j  ’ 


.HE  CHE 


hHE  fHE  H(t) 


cot  $  . 


When  x  ^_0.25,  the  HMD  expressions  listed  in  the  previous  section  reduce 


mn„G(t)  cos  Pe_ap 
PHM  _ 


mn_H(tj  sin  pe' 


i^u.mG(u)  cos  pe" 
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uHM  mH(t)  sin  4>e“ac 

"°  2^ 

UHM  mG(t)  cos  $e~ 


in0kpG(u) 

G(t) 

inQkpG(u) 
neG(t)  ’ 


hHM  PHM 

_1_  .  .  G(tl  cot 

hHM  eHM  H(t) 

P  $ 

When  p  <  h/3,  G(u)  -  G(t)  -  V(t)  -  1.0,  H(t)  -  2,  and  ap  -  0.  For  this 
case,  the  formulas  presented  in  this  section  [equations  (56)  through  (76)), 
reduce  to  the  familiar  quasi -nearfield  range  results. 5-8  They  will  not  be 
repeated  here  since  they  are  already  given  by  equations  (56)  through  (76)  with 
G(u)  -  G(t)  -  V(t)  -  1.0,  H(t)  -  2.0,  and  ap  =  0. 


FIELD-STRENGTH  FORMULAS  FOR  c  >  5h 


When  p  >  5h, 


/->  f  \  wU  U 

G(u) - =  -  , 

TT  h 


G (t)  -  H(t)  -  — 

ir 


h(c/v)' 


For  this  case,  the  general  field-strength  formulas  presented  in  equations 
(31)  through  (53)  reduce  to  those  of  Wait*  and  Galejs.-  For  the  sake  of  com¬ 
pleteness,  they  will  be  repeated  here. 
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FIELD-STRENGTH  FORMULAS  FOR  p  >  3h  AND  x  >  1.6 


For  x  >  1.6,  the  Hankel  functions  can  be  approximated  by 


H0(z) (x)  -  . 

(101) 

Hp>W  -  i^e-i(X-*/4)  , 

(102) 

(■l")HlU)t*)  ->/Ye-i<X"/4)  . 

(103) 

(•"jKp’wfM  -  ix^e-K*-*/4)  . 

(104) 

Thus,  the  VED  expressions  listed  in  the  previous  section 

reduce  to 

-otp 

VE  nePe  Fx  ifx-ir/4) 

be  '  '  2irhp  \  2 

(105) 

VE  noP(c/v)e-™  _  x.,/4) 

-  '  ’  2whp  \  2  6 

(106) 

H;E  •  PC^e'1(X‘,/J)  • 

HO") 

eVe 

^E  '  -  n0(G/v)  ’ 

(108) 

eve 

~  -  -  (n0/ne)  (c/v)  . 

E 

(109) 

'.'.hen  x  l.t>,  the  HED  expressions  listed  in  the  previous 

section  reduce 

j-HE  P  cos  $e’ao  ( ix)  Fx  Q-i  (x-tt/4) 

2*t  ( oe  +  iuj£e)h  (c/v)  2o2  \  2 

,-HE  p  sm  :e'10  Fx  -iix-~/4  ! 

2-(-e  +  i  -  e  )  h  (c/v)-w-\- 

(110) 

(111) 

ail) 


TR  7521 


ilOUnP  cos  Ae_otp  , — 

_  [HE  e-i(x-ir/4) 

2iryehp  V  2 
p  sin  4>e'ap  px  -i(x-ir/4) 

2iTYeh(c/v)2p2,V  2 


hHE  _  _  P  cos  4>e~ap  (ix)  fnx  -i  Cx-ir/4) 
^  2TrYeh  (c/v)  2p  2  \  2 


-n0 (c/v)  , 


EHE 

Je  '  Cn0/ne3  Cc/v)  , 

HP 


and 


hhe 


ix  cot  <J) 


[sin(o/a)  I 
(P/a)  J  ‘ 


(11 

(11 

(11. 

(11 

(lli 

(11 


to 


When  x  >  1.6,  the  HMD  expressions  listed  in  the  previous  section  reduce 


mn„  cos  ie'ap(ix)  , 
eH*'«  .  e  p*  e-i(x--/4) 


2-h(c/v)-- 


(lli 


-mn  sxn  :e  , 

:HM  .  JL _  jrx  -“i  (x-t/4  J 


J  2TTh(c/v)2o2  *  2 
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and 


hhm 


ix  cot  $ 


[sin(p/a) 

(P/a) 


]• 


(125) 


DISCUSSION 


We  have  used^  the  recently  developed  theory  of  Greifinger  and  Greifin- 
ger10-12  anci  the  Wait  exponential  ionospheric-conductivity  profile  to  deter¬ 
mine  average  ELF  propagation  constants  for  both  daytime  and  nighttime 
propagation  conditions.  The  resulting  average  values  of  ELF  attenuation  rate, 
phase-velocity  ratio,  and  ionospheric  reflection  height  (presented  in  table  1) 
are  in  excellent  agreement  with  measured  data.® 

To  shed  further  light  on  the  nature  of  the  ELF  field  strengths  in  the 
earth-ionosphere  waveguide,  the  radial  impedance  of  the  wave  (Ej/H^)  now  is 


Table  1.  Typical  ELF  Propagation  Parameters 


(• 
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considered.  By  definition,  this  quantity  is  the  impedance  of  the  wave  looking 
in  the  radial,  or  p,  direction. 

For  x  >>  1 , 

E, 


Hi  '  -notc/v)  ’ 


while,  for  x  <<  1  and  o  <  h/3. 


.VE  ikp 


gHE  gHM 

TmE  '  '  -lk2n0 

Hj  H> 


For  the  intermediate,  and  most  interesting,  range, 


in0(c/v) 

V(t)  ♦  G(u)x2hJ2)(x) 

xG(u) 

^■h(2)(x) 

EfJE 

ehm 

i " - xG ( u ) 

5? 

hHm 

(c/vjf(x)G(.t  j 

while,  for  x  <_  0.25,  equations  (33)  and  (42)  reduce  to 
cVE 


'2  M0 
>VE  (ko)G(u) 


{vet) 


c(u)x;[ln(iJL5i)  .  i2]j 


,-HL  .  ilM 

E  -  i. . 

J4E  '  JiM 
H;  hj 


1  i\ ..  kj  i  U 

G(t) 


(126) 


(127) 


(123) 


Referring  to  equation  (59),  we  see  that  the  minimum  value  of  e\E/H^E 
occur  when  V  t  ,  G • u . x  -ns  1 . 123  \ ; .  That  is , 


V  t- 


r.' 

r»  .* 
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ti 
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cVE 


jVE 


(120tt)  (rr/2)  (kp)  (c/v): 


MIN 


=  60ir2  (kp)  (c/v)  2  . 


(129) 


Alternatively,  the  approximate  distance  from  the  VED  source  where  the 
minimum  value  of  radial  wave  impedance  occurs,  omv,  can  be  expressed  as  (from 
equation  (129)) 


80.63  E^/H^ 

z  4> 


nun 


"mv 


f (c/v)2 


km  . 


(130) 


Presented  in  figures  5  and  6  are  plots  of  the  VED  radial  wave  impedance 
versus  distance  for  frequencies  of  30  to  300  Hz.  Equation  (33)  and  the  values 
of  h  and  c/v  listed  in  table  1  were  used  in  the  calculations .  Note  that,  for 
frequencies  of  30  to  100  Hz,  there  is  a  unique  distance  where  the  minimum  value 
of  radial  wave  impedance  occurs.  Presented  in  table  2  are  values  of  c  cal¬ 
culated  from  equation  (130) .  Comparing  these  values  with  the  curves  of  figures 
5  and  6  reveals  that  the  table-2  omv  calculations  are  accurate  within  10  km. 

Presented  in  figures  7  and  8  are  plots  of  the  HED  and  HMD  radial  wave 
impedance  versus  distance  for  frequencies  of  30  to  300  Hz.  Equation  (42)  and 
the  values  of  h  and  c/v  listed  in  table  1  were  used  in  the  calculations. 


Table  2.  Approximate  Distance  Where  the  Minimum  Value 
of  the  VED  Radial  Wave  Impedance  Occurs 
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Referring  to  figures  5  through  8,  we  see  that  there  is  a  substantial  var¬ 
iation  in  both  the  VED  and  HED  radial  wave  impedance.  For  example,  at  30  Hz, 
the  VED  wave  impedance  is  equal  to  30,000  ohms  to  20  km,  120  ohms  at  180  km, 

and  505  ohms  (120irc/v)  at  5,000  km  (figure  5).  On  the  other  hand,  the  30-Hz 

HED  (or  HMD)  wave  impedance  varies  from  5  ohms  at  20  km  to  505  ohms  at  5,000 

km  (figure  7) . 

As  Waitl  has  pointed  out,  the  observed  variation  of  the  magnitude  and/or 
phase  of  the  radial  wave  impedance  as  a  function  of  frequency  should  provide  a 
basis  for  distance  measuring  (provided  the  atmospheric  source  could  be  repre¬ 
sented  by  an  equivalent  VED  or  HED).  Such  a  scheme,  while  admittedly  crude, 
requires  only  one  receiving  station  equipped  with  a  vertical  whip  and  loop 
antenna. 


CONCLUSIONS 


In  this  report,  we  have  presented  new  formulas  for  HED,  HMD,  and  VED  ELF 
radio-wave  propagation  in  the  earth-ionosphere  waveguide.  These  new  formulas 
extend  the  results  of  Wait^  and  Galejs,-  which  are  valid  for  measurement  dis¬ 
tances  greater  than  approximately  three  ionospheric  reflecting  heights,  down 
to  the  quasi-nearfield  range,  which  is  defined  as  the  range  where  the  measure¬ 
ment  distance  is  greater  than  an  earth  wavelength  but  much  less  than  a  free- 
space  wavelength.  For  the  sake  of  completeness,  the  abovementioned  previously- 
derived  formulas  also  have  been  included. 

Plots  of  the  VED,  HED,  and  HMD  radial  wave  impedance  versus  distance  have 
been  presented  for  both  daytime  and  nighttime  propagation  conditions.  These 
plots  show  a  substantial  variation  in  the  radial  wave  impedances  for  distances 
less  than  1,100  km.  Also,  we  have  shown  that  there  is  a  unique  distance  where 
the  minimum  value  of  the  VED  radial  wave  impedance  occurs . 


|Ez/H*|  (ii) 


10 


100 


DISTANCE  a  (km) 


1 ,000 


Figure  7.  HED  and  HMD  Radial  Wave  Impedance  Versus  Distance 
for  Daytime  Propagation  Conditions 


DISTANCE  a  (km) 


Figure  S .  ilLL>  and  HMD  Radial  Wave  Impedance  Versus  Distance 
for  Nighttime  Propagation  Conditions 
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Appendix 

E/H  CALCULATION  PROGRAMS 

Five  programs  were  written  to  solve  for  E/H  fields  at  various  distances 
from  a  transmission  source.  The  following  cases  were  considered: 

1.  Ez/H-p  of  an  HED  source,  nearfield  approximation; 

2.  Ez/H^  of  a  VED  source,  nearfield  approximation; 

3.  Ez/Hp  of  an  HED  source,  nearfield  approximation; 

4.  Ez/H^  of  an  HED  source,  farfield  approximation;  and 

5.  E,/Hj  of  a  VED  source,  farfield  approximation. 

The  programs  are  written  in  VAX-11  FORTRAN.  They  reside  as  executable 
modules  on  the  Naval  Underwater  Systems  Center  (NUSC)  VAX  system.  They  can  be 
invoked  by  keying  RUN  'filename'  where  the  five  files  shown  below  correspond 
to  the  five  cases  listed  above: 

1.  v703: :DU01: [ETC.EHNEARFLDJEHHED; 

2.  v703: : DUO 1 : [ ETC . EHNEARF LD ] EHVED ; 

3.  v“05: : DU01 : [ETC . EHNEARF LD ] E2HRH0 ; 

4.  v703 : : DU01 : [ETC . EHFARF LD J  EHHED ;  and 

5.  v703: :DU01: [ETC. EHFARFLD] EHVED. 

Each  of  the  programs  contains  an  interactive  section  in  which  the  user 
enters  parameters  of  interest.  These  parameters  involve  the  range  of  distance 
from  the  source  for  which  E/H  is  to  be  calculated,  the  frequency  of  interest, 
the  time  of  day,  and,  sometimes,  the  angle  from  the  cenith. 

Listings  are  provided  for  each  of  the  programs  created.  Program  structure 
and  flow  are  clearly  indicated  in  the  listings.  Most  of  the  program  varia¬ 
bles,  such  as  T,  C.T  i=  G ,  1 1  .  ,  X,  R1I0,  Rl'I.HT  (=  ionospheric  reflection  height', 
CYEE  ,  =  c  v  ,  TOD  ■ =  night  or  day  values,,  etc.,  correspond  exact ly  to  the 

variables  in  the  equations.  It  is  suggested  that  the  reader  have  the  equa¬ 
tions  at  hand,  particularly  when  tracing  the  flow  of  the  actual  calculations , 
as  the  calculations  follow  the  equations  in  a  logical  and  consistent  manner. 
These  programs  are  pretty  much  'number  crunchers;'  their  flow  is  not  complex. 

Calls  are  made  for  calculation  of  hyperbolic  functions  and  Hankel  func¬ 
tions.  Hie  hyperbolic  function  routines  also  are  written  in  \A.\-11  FORTRAN; 
the  listings  for  these  routines  have  been  included.  The  Hankel - funct ion 
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calculation  consists  of  calling  the  appropriate  IMSL  Bessel  functions  and  com¬ 
bining  them  in  the  calling  program,  as  the  Hankel  function  is  a  combination  of 
Bessel  functions.  The  IMSL  routines  are  available  through  the  'IMSLIBS' 
library  on  the  VAX. 


CREATED: 

LAST  UPOATE; 
BY: 

PURPOSE: 


27FEB85 
07MAR85 
A.  KUZEL 

FORTRAN  PROGRAM  TO  CALCULATE  THE  VALUE 
OF  E/H  FOR  AN  ELF  WAVEFORM  OF  A  SET  OF 
ALLOWED  FREQUENCIES  F  AT  VARIOUS 
DISTANCES,  <■  200  KM.  FROM  A  HE D  SOURCE. 
THIS  IS  THE  NEAR  FIELD  APPROXIMATION. 


•OECLARE  VARIABLES* 


COMPLEX 

REAL 

real 

CHARACTER 


I , HQ , H 1 , HR , RAZ 

K , JO . J 1 , P I . RHO . FREQ . I NC . MAXDI S . T , U . V , X , GT , VO ( 1 ) 
CVEE , RFLHT , MMBS JO . MMBSJ 1 , RAZA , RAZPH , Y 1 ( 1 ) 

TOO 


•INITIALIZE  constants* 


K=3.00E-05 
PI=4*ATAN( 1.0) 

I -CMPLX (0.0, 1.0) 


THE  FOLLOWING  SECTION  IS  THE  INTERACTIVE  PORTION  OF  THE  PROGRAM. 
THE  USER  CHOOSES  STARTING  DISTANCE.  INCREMENT.  AND  ENDING  DISTANCE 
TO  CALCULATE,  NIGHT  OR  DAY  VALUES.  AND  THE  FREQUENCY  OF  INTEREST. 


WR I TE ( 5 , 200 ) 

READ ( 5 . 204 ) RHO 
WR I TE l 5 , 20  1  ) 

READ ( 5 , 204) INC 
WRITE(5 . 202) 
READ(5.204)MAXDIS 
WR I TE ( 5 , 203 ) 
REAO(5,204)FHEQ 
WR I TE ( S , 205  ) 
R£AD(5,206)TOD 


•SELECT  PROPER  IONOSPHERIC  REFLECTION  HEIGHT  AND 
C/V  CONSTANTS  BASED  ON  USER'S  INPUT  DATA* 


IF(FRE0. EQ. 30)  THEN 
I F ( TOO . EQ . 'O' )  THEN 
RFLHT=46 . 1 
CVEE*  1 . 34 
END  IF 

I F ( TOO . EQ . 'N' )  THEN 
RFLHT  =72. 0 
CVEE* 1.12 
ENO  IF 
ENO  IF 

I F ( FREQ . EQ . 50 )  THEN 
I  F  (  TOO  .  EQ  .  0  )  THEN 
RFLHT*47 . 8 
CVEE*  1 .30 
ENO  IF 

I F ( TOO . EQ . ' N ’ )  THEN 


V.V.s.V, 

w.  *  J * 
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I F (TOO . EQ . ' N ' )  toft  I  TE  ( 10,212) 
toft ITE ( 10 , 20B ) 
toft ITEC 10 , 209 ) 

WRITE( 10,208) 

•PERFORM  THE  ACTUAL  CALCULATION* 

X=2»PI*FREQ»RH0«CVEE/K 

T*RH0*PI/(2»RFLHT»CVEE»»2) 

U*T 

V»T 

CALL  COTH(U) 

CALL  CSCH(V) 

GT»2«T/PI*U 

GT*GT*-(  1-2/PI )  *T**2*V**2 
JO«MMBSJO(X. IER) 

J  1*MM8S J 1 (X , IER) 

CALL  MMBSVN ( X.O.  ,  1 , VO , I ER ) 

CALL  MMBSVN (X.O. 999, 1. VI, IER) 
HO=CMPLX<  JO ,-V0( 1) ) 
H1*CMPLX(J1.-V1(  1)) 

HR=M0/M1 

RAZ=l-X»Hft 

RAZ=1/RAZ 

RAZ3RAZ*120*PI*I*X 
RAZ=-RAZ 
U=RFLHT/RhO 
CALL  TANM(U) 

RAZ*RAZ/(CVEE»GT*U) 

RAZA=CABS(RAZ) 

RAZPH*ATAN2(AIMAG(RAZ) .REAL(RAZ) ) 

RAZPM*RAZPH* 180 .00/PI 

toR  I  TE ( 10, 210)RMO,RAZa,RAZPH 

RHO=RHO*INC 

IF(RmO.LE.MAXDIS)  GO  TO  400 
•FORMAT  STATEMENTS* 


FORMAT! 4X , 
FORMAT ( J* , 
FORMAT ( 4X  . 
FORMAT ( 4X , 
FORMAT ( F 1 2 
FORMAT ( 4X , 
FORMAT ( A  2 ) 
FORMA  T ( 9X , 
FORMAT ( 5X , 
FORMATt 1 IX 
FORMAT ( IOX 
FORMAT ( 9X , 
FORMA  T ( 9X . 
FORMAT ( 9x . 
STOP 
ENO 


ENTER  STARTING  DISTANCE  FROM  SOURCE;  ) 
ENTER  DISTANCE  INCREMENT;) 

ENTER  MAXIMUM  DISTANCE  TO  COMPUTE;’) 
ENTER  FREQUENCY : ' ) 


NIGHT  (N)  OR  DAV  ( D ) 7 ; ' ) 

HORIZONTAL  electric  dipole 


FREQ  =  ' ,F7.2) 


’OISTANCE’ , 8X . -MAGNITUDE’ . 9X , ’PHASE’ ) 
F 7 . 2 ,  IOX, F7 .  2 ,  IOX , F7 . 2) 

TOO  =  DAYTIME ’ ) 

TOO  =  NIGHT’ ) 

E/H  NEAR  FIELD  APPROXIMATION') 


n  n  n  n  n  r>  o  n  o  n  o  o  o  o  o  r»  o 
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CREATED :  28FEB85 

LAST  UPOATE:  07MAR85 

BY:  A.  KUZEL 

PURPOSE:  FORTRAN  PROGRAM  TO  CALCULATE  THE  VALUE 

OF  E/H  FOR  AN  ELF  WAVEFORM  OF  A  SET  OF 
ALLOWED  FREQUENCIES  F  AT  VARIOUS 
DISTANCES.  <=  200  KM.  FROM  THE  VED  SOURCE. 
THIS  IS  THE  NEAR  FIELD  APPROXIMATION. 


•OECLARE  VARIABLES* 

COMPLEX  I .HO.Hl.HR.RAZ 

REAL  K. JO. Jl.PI , RHO , FREQ , I NC , MAXDI S , T , U , V , X , GT , VO ( 1 ) 

REAL  CVEE . RFLHT , MMBS JO . MMBS J 1 , RAZA . RAZPH . Y 1 ( 1 ) 

CHARACTER  TOD 

•INITIALIZE  CONSTANTS* 

k  =  3 . OOE  *05 
P I =4* AT  AN ( 1.0) 

I-CMPlX(0.0. 1.0) 


The  FOLLOWING  SECTION  IS  THE  INTERACTIVE  portion  OF  THE  PROGRAM  • 
The  USER  CHOOSES  STARTING  DISTANCE.  INCREMENT.  AND  ENOING  DISTANCE  • 
TO  CALCULATE.  NIGHT  OR  OAV  VALUES.  AND  FREQUENCY  OF  INTEREST.  • 


WR I TE ( 5 , 200 ) 

REA0(5. 204JRHO 
WR I TE ( 5 . 20 1 ) 

REAO(5 . 204) INC 
WR I TE ( 5 . 202 ) 

REAO( 5 . 204)MAXOI S 
WR I TE ( 5 , 203 ) 

R  E A  0 ( 5 . 204 ) FREQ 
*RI T£ ( 5 . 205 ) 

REAOl 5 . 206 i '00 

C  *CHOOSE  PROPER  VALUES  OF  IONOSPHERIC  REFLECTION  HEIGHT  AND 

C  C/V  CONSTANTS  BASED  ON  USERS  INPUT  DATA* 


IF (FREQ . EQ . 30 )  THEN 
I F ( TOD . EQ .  O’ )  THEN 
RFLHT=46 . 1 
CVEE  *  1 . 34 
ENO  IF 

I F ( TOO . EQ . ' N ’ )  THEN 
RF  LHT  =  72. 0 
CV£E=  1  .  12 
ENO  r  F 
ENO  IF 

I F ( FREQ . EQ  50 )  THEN 
I F ( *00 . EQ .  0  i  ’HEN 

RFlmT-47  8 
CvEE= 1 . 30 

ENO  IF 

IF (TOO  EQ.  N  )  THEN 
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RFLHT=73 . 3 
CVEE* 1 . It 
END  IF 
END  IF 

I F ( FREQ . EQ . 75 )  THEN 
I F (TOD . EQ . 'O' )  THEN 
RFLHT  =  49 . 1 
CVEE= 1.27 
END  IF 

I F ( too . EQ . ' n ' )  then 
RFLHT=74.3 
CV£E»1. 10 
END  IF 
END  IF 

IF(FREQ.EO. 100)  THEN 
IF(TOO.EQ. '0' )  THEN 
RFLHT =50 . 1 
CVE£=1 . 25 
ENO  IF 

I F (TOD . EQ . 'N' )  THEN 
RFLHT=7S . 0 
CvE£= I . 09 
ENO  IF 
ENO  IF 

IFCFREO. EQ. 150)  THEN 
I F ( TOO . EQ . 'O' )  THEN 
RFLHT  =5 1 . 4 
C VEE= 1.22 
ENO  IF 

I F ( TOO . EO . 'N' )  THEN 
RFLHT=76 . 0 
CVEE= 1 . 09 
ENO  IF 
ENO  IF 

IF(FREO.EQ.200)  THEN 
I F ( TOO . EQ .  O' )  Then 
RFlhT  =  52 . 4 
CVEE=  I  . 20 
END  IF 

;F(TOO.£0.  N  )  'MEN 
RFlhT  =  76 . a 
CVEE=1.08 
ENO  IF 
ENO  IF 

IF(FRE0. EQ. 300)  Then 
I F ( TOO . EQ .  O' )  THEN 
RFLHT  =  53 . 7 
CvEE  =  1  .  18 
ENO  IF 

IF(TOO.EQ.  N' )  THEN 
RFlhT  =  77 . 8 
CVEE*  1.07 
ENO  IF 
ENO  IF 

•WRITE  HEADERS  FOR  OuTPuT  TABLE* 

wR I TE i  10.220  ) 

WRITE! 10 . 208 ) 

WRITE!  10 , 207 ) FREQ 
I F ( TOD . EQ .’O')  WRITE! 10.211) 


1  F  (TOO  .  EQ  .  '  N  '  )  MR  I  TE (10,212) 

MR ITEC 10.208) 

MR  I TE ( 10 , 209 ) 

MR  I TE ( 10 , 208 ) 

C  ‘PERFORM  THE  ACTUAL  CALCULATION* 

400  X=2»PI •FREQ»RHO»CVEE/K 

T  =  RHO*PI/  (  2»RFLHT*CVEE“2) 

U=T 
V  =  T 

CALL  COTH(U) 

CALL  CSCH(V) 

GT=2»T/PI »U 

GT=GT*(  1-2/PI  )»T“2»V»»2 
VT=T“3*V**2*U 
HT*VT*GT 
J0=MMBS JO ( X , IER) 

J 1=MMBS J 1  ( X , IER) 

CALL  MMBSVN(X .0 . . 1 . VO. IER) 

CALL  MMBSVN ( X. 0.999, 1, VI, IER) 

HO  =  CMPLX( JO . -VO (  1 ) ) 

H1=CMPLX( j 1 , -V 1(  1 ) ) 

RA2=h0*X“2*CVEE“2‘HT‘  (  PI  /  2  )  *1 
RA2=RA2*VT 

RA2=RA2/ (Hl«X»Pt • ( -I ) / 2 ) 

U=RFLHT/RHO 
CALL  TANH(U) 

RA2-RA2*U/X 

RA2fRA2‘I*120*PI‘CV££ 

RA2A=CABS(RA2) 

RA2PH=ATAN2(AIMAG(RA2) ,REAL(RA2) ) 

RA2Ph=RA2PH* 180 . 00/PI 

MR  I TE ( 10 , 210) RHO , RA2a , RA2PM 

RH0=RH0*INC 

if(rho.le.maxdis)  go  to  400 
•format  STATEMENTS* 

FJRMAT(4x,  ENTER  STARTING  DISTANCE  FRGM  SOURCE:! 

201  F0RMAT(4X.  ENTER  OISTANCE  INCREMENT:) 

202  FORMAT (4X, 'ENTER  MAXIMUM  DISTANCE  TO  COMPUTE:) 

203  FORMAT (4X, -ENTER  FREQUENCY : ' ) 

204  FORMATfF 12 . 4) 

205  F0RMAT(4X, 'NIGHT  (N)  OR  DAY  ( D ) 7 : ' ) 

206  FORMAT ( A2 ) 

207  FORMAT ( 10X, 'VERTICAL  ELECTRIC  OIPOLE  FREQ  =  -,F7 

208  F0RMAT(5X, •  ’ ) 

209  FORMAT ( 1 IX ,  DISTANCE' ,8X. 'MAGNITUDE' , 10X, ’PHASE' ) 

210  FORMAT ( 10X , F 7 . 2 , 9X , F9 . 2 , 10X , F 7 . 2 ) 

211  FORMAT ( 1 2X , ' TOO  =  DAYTIME') 

212  FORMAT! 12X ,- TOD  =  NIGHT') 

220  F0RMAT(9X, "E/H  NEAR  FIELD  APPROXIMATION’) 

STOP 

END 


o  n  o  o  o  r>  o  o  oonooooonno 
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CREATED: 

LAST  UPDATE: 
BV: 

PURPOSE: 


01MARSS 
01MAR85 
A.  KUZEL 

fortran  program  TO  CALCULATE  The  value 
OF  E  SUB  Z  OVER  H  SUB  RhO  FOR  AN  ELF 
WAVEFORM  OF  A  SET  OF  ALLOWED  FREQUENCIES  F 
AT  VARIOUS  DISTANCES.  <  =  200  KM.  FROM  AN 
HED  SOURCE.  AT  VARIUOS  ANGLES  FROM  THE  ZENITH. 
THIS  IS  THE  NEAR  FIELD  APPROXIMATION 


•DECLARE  VARIABLES* 


COMPLEX 

REAL 

REAL 

CHARACTER 


I  .  RAZ 

K.PI .RHO, FREQ. I NC . MAXDI S . T , U . V , X , GT 
CVEE .RFLHT, RAZ A , RAZPH.VT.HT. PHI 
TOO 


•initialize  variables* 


K  =  3 . 0OE»O5 
PI =  4* ATAN (1.0) 
I=CMPLX(0 .0 . 1.0) 


THE  FOLLOWING  SECTION  IS  the  INTERACTIVE  PORTION  OF  THE  PROGRAM. 
THE  USER  CHOOSES  STARTING  OISTANCE.  INCREMENT.  AND  ENDING  DISTANCE 
TO  CALCULATE.  NIGHT  OR  OAV  VALUES.  FREQUENCY  OF  INTEREST.  AND  ANGLE 
OF  INTEREST. 


WRITE(5. 200) 

READ ( 5 . 204 ) RHO 
WR I TE ( 5 . 20  1  ) 

REA0(5,204)InC 
wR  I  TE  (  5 , 202  ) 

REAQ(5, 204 )MAaDI S 
.  *£»  ’<•**  5.203) 

READ: 5 . 204 ; FREQ 
WRI TE(5 . 205  ) 

REAO(5.206)TOO 
WRITE(5 . 213 ) 

REAOC5. 204 ) PHI 

C  ‘SELECT  PROPER  IONOSPHERIC  REFLECTION  HEIGHT  AND 

C  C/V  VALUES  BASED  ON  USER'S  INPUT  DATA.* 

I F  < FREQ . EQ . 30 )  THEN 
I F ( TOO . EQ . ' D ' )  THEN 
RFLHT  =  46 .  1 
C VE  E  =  1  . 34 
ENO  I  F 

IF ( TOD .  EQ  . ' N  )  THEN 
RF  lht  =  7  2 .0 
Cv  £  E -  1  12 

ENO  I  F 
END  I  F 

I  F  (  FREQ  .  EO  .  50  )  Tr-EN 
I F ( TOD . EQ .  0  )  THLN 
RFLHT  =  47 . 8 
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cvee=i.30 
ENO  IF 

I F (TOD . EQ . 'N- 
RFLHT=73.3 
CVEE=l. 11 
ENO  IF 
END  IF 

I  F ( FREQ . EQ . 75) 

I F ( TOO . EQ . 'O' 
RFI_HT  =  49 . 1 
CV££= 1 . 27 
ENO  IF 

I F (TOO . EQ . 'N' 
RFLHT=74 . 3 
CVEES 1 . 10 
ENO  IF 
ENO  IF 

IF(FREO.EQ.  100) 
I F (TOO . EQ . 'O' 
RFLMT=SO . 1 
CvEE  =  I . 25 
ENO  IF 

IF (TOD. EQ. 'N ' 
R  F  LH  T  -  7  5 . 0 
CvE£= 1 .09 
END  IF 
ENO  IF 

I F ( FREQ . EQ . ISO) 
I F ( TOO . EQ . 'O' 
RFLHT=51.4 
CVEE* 1 . 22 
END  IF 

I F ( TOO . EQ .  N- 
RFU1T*76.0 
cvee= i,09 

ENO  IF 
ENO  IF 

IFI  FREQ  .  EQ.  200  ■) 
IF.  TOO.zQ.  0 


)  THEN 


THEN 

)  then 


)  then 


then 

)  THEN 


)  Then 


then 
)  then 


)  THEN 


thEn 

)  'HEN 


2 vE E  =  l  . 20 
ENO  IF 

IF(TOO.£Q.  N' )  THEN 
RF  LHT  =  76 . 8 
CVEE= 1 . 08 
END  IF 
ENO  IF 

IF( FREQ . EQ . 300 )  THEN 
IF(TOO.£Q. 'D  )  THEN 
RFi_HT  =  53 . 7 
C VEE 1 1  .  18 
ENO  IF 

I  F ( TOO . EQ .  N  ' )  Then 
RFLMT -77 . 8 
CvEE  = 1 .07 
END  IF 
END  IF 


•wRITE  headers  for  Juipul  facSlE* 


wR I TE (  10 , 207 )FR£Q 
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I F (TOO . EQ .'O')  MR  I TE ( 10,211) PHI 
IF (TOD. EQ.  ' N ' )  MR  I  TEC  10, 2 12} PHI 
MR I  TEC  10,206) 

MR  I TE  C 10 , 209 ) 

MR  I  TEC  10,206) 

PHI=PHI‘PI/180.00 

C  ‘PERFORM  THE  ACTUAL  CALCULATION* 

400  X=2‘PI‘FREQ»RH0‘Cv£E/K 

T =RHO‘PI / C  2‘RFLHT  »CVEE“ 2 ) 

U*T 

V*T 

CALL  COTHCU) 

CALL  CSCHCV) 

GT=2‘T / PI  *U 

GT=GT*C  l-2/PI)*T“2‘V“2 

VT=T“3‘V“2‘U 

HT=vT»GT 

U=( 1/TANCPHI ) ) 

V=RFlhT/RhO 
CALL  TANH(V) 

RA2=I • 12Q‘PI »X‘U 
RAZ=RAZ/(CvEE*hT‘V) 

RAZA-CABSCRA2) 

RA2PH=AT  AN2 ( A  I MAGC  RA2 ) .REALCRAZ) ) 

RAZPH=RAZPH‘ 180 . 00/PI 

MR  I  TEC  10,210) RHO . RA2A , RAZPH 

RHO=RHO*INC 

IFCRHO.LE.MAXOIS)  GO  TO  400 
C  ‘FORMAT  STATEMENTS* 

200  F0RMATC4X. ’ENTER  STARTING  DISTANCE  FROM  SOURCE:') 

201  F0RMAT(4X, ' ENTER  DISTANCE  INCREMENT:') 

202  F0RMAT(4X , • ENTER  MAXIMUM  DISTANCE  TO  COMPUTE:) 

203  F0RMAT(4X, -ENTER  FREQUENCY:') 

204  FORMAT ( F 12 . 4 ) 

225  F0RMATC4X, 'NIGHT  (N)  OR  DAY  ( 0 ) ? :  ) 

2  jo  format ( A2 ) 

207  FORMAT ( 10X , ■ MED  SOURCE  EZ/HRhO  FREQ  =  '.F7.2) 

208  FORMAT (SX, '  ' ) 

209  FQRMATC 11X, ’DISTANCE' ,8X, 'MAGNITUDE' , 10X, 'PHASE' ) 

2 10  FORMAT ( 10X , F7 . 2 ,9X, F9 . 2 . 10X.F7 . 2 ) 

211  FORMATC 12X , ' TOO  =  OAYTIME  PHI  =  '.F5.2) 

212  FORMATC 12X. 'TOO  =  NIGHT  PHI  =  ' , F5 . 2 ) 

2 1 3  F0RMATC4X, 'ENTER  PHI  IN  DEGREES:') 

STOP 

END 


\-  1  ■* 


oooononooo 
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OATE  CREATED:  2 1FEB8S 
LAST  UPOATE:  07MARBS 

BV:  A.  KUZEL 

PURPOSE:  THIS  PROGRAM  CALCULATES  THE  VALUE  OF  THE  E/H 

FIELD  OF  AN  ELF  WAVEFORM  OF  ONE  OF  A  SET  OF 
FREQUENCIES  AT  VARIOUS  DISTANCES  >=  200  KM 
FROM  AN  HED  SOURCE.  THIS  IS  THE  FAR  FIELD 
APPROXIMATION. 


C  ‘DECLARE  VARIABLES* 

COMPLEX  I. V. HO. HI, HR 

REAL  JO. Jl.VOC 1) .Yi( 1) , MMBS JO , MMBS J 1 . MAXDI S 

REAL  RHO.INC.PI , CVEE . FREQ , FUDGE . X , VA . PH 

CHARACTER  TOD 

c  ‘Initialize  constants* 

PI=4‘ATAN(  1.0) 

I =CMPLX (0.0, 1.0) 

. . . . . . . 

C  the  FOLLOWING  SECTION  IS  The  INTERACTIVE  PORTION  OF  THE  PROGRAM 

c  the  user  chooses  starting  oistance.  increment,  and  ending  distance 
c  TO  calculate,  night  or  oav  values,  AND  The  frequencv  of  interest. 

C  “  . “ . . . . 

WR I TE ( S . 540 ) 

REA0(5.545)RH0 

WRITE(S,547) 

READ( 5 , 560 ) I NC 
FuDGE=3.0E-05 
wR I TE ( 5 , 548 ) 

READ( 5 . 545 )MAX0I S 
wR I TE ( 5 . 550 ) 

R£A0( 5 , 560) FREQ 
wRI TE( 5 , 205  I 
R  E A 0 ( 5 , 206  >  T  30 

C  ‘SELECT  PROPER  IONOSPHERIC  REFLECTION  HEIGHT  AND 

C  C/V  CONSTANTS  BASEO  ON  USERS  INPUT  DATA* 

I F ( FREQ . EQ . 30 )  THEN 
I F ( TOO . EQ . ’O' )  THEN 
RFlhT  =  46 .  1 
CVEE= 1 , 34 
ENO  IF 

I  F  (  TOO  .  EQ  .  -  N  '  )  THEN 
RFlhT  =  72 . 0 
CVEE  = 1.12 
ENO  IF 
ENO  IF 

IF(FREQ . EQ . 50 )  THEN 
I F ( TOO . EQ .  D  )  Then 
RFlhT  =  47  .  8 
C vEE  = 1  30 
ENO  IF 

I F ( TOD . EQ . ’N' )  THEN 
RFLHT  =  73 , 3 
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CV££=1. 11 
END  IF 
ENO  IF 

I F (FREQ . EQ . 75 )  THEN 
I F(TOD . EQ . 'O' )  THEN 
RFLHT =49 . 1 
CV£E  = 1.27 
END  IF 

I F (TOO . EQ . ‘N- )  Then 
RFLHT=74. 3 
CVEE=1 . 10 
ENO  IF 
ENO  IF 

IF (FREQ. EO. 100)  THEN 

I f ( too . EQ . '0-)  then 

RFLHT =50 . 1 
CVEE* 1 . 25 
ENO  IF 

I F ( TOD . EQ . 'N- )  THEN 
RFLHT  =  75.0 
Cv£E= 1 . 09 
ENO  IF 
ENO  IF 

IF(FREQ.£Q. 150)  THEN 
I F ( TOO . EQ . '0' )  THEN 
RFLHT -5 1 . 4 
CVEE- 1.22 
ENO  IF 

I F ( TOO . EQ . ' N ' )  THEN 
RFLHT=76.0 
CVEE= 1 . 09 
ENO  IF 
ENO  IF 

IFCFREQ.EQ.200)  THEN 
I  F  (  TOO  .  EQ  .  ’0  )  THEN 
RFLHT  =  52 . 4 
CVEE=  1  . 20 
ENO  IF 

I F ( TOO . EQ . ' N ' }  THEN 
RFlhT  =  76 . 8 
CVEE  = 1 . 08 
ENO  IF 
ENO  IF 

IF(FREQ. EQ. 300)  THEN 
I F ( TOD . EQ . '0)  THEN 
RFlhT  =  53 . 7 
CVEE=1. 18 
END  IF 

I F (TOO . EQ .  N' )  THEN 
RFLHT  =  77 . 8 
CvEE  = 1.07 
ENO  IF 
ENO  IF 

JF(FR£Q.EQ.400)  CVEE  =  1.15 
I F ( FREQ . EQ . 800 )  CVEE=1.  11 
I F ( FREQ  EQ . 1600)  CVEE=1.07 

•write  heaqers  for  output  table* 

WRITE( 10,570)FREQ 
IFCTOD.EQ. 'O' )  WR ITEC10.211) 
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I F (TOO . EQ .  ' N ' )  MR  I TE ( 10 , 2 12 ) 

MR  I TE ( 10 , 575 ) 

MR  t TE ( 10 , 580 ) 

WRITE( 10.575) 

•PERFORM  THE  ACTUAL  CALCULATION* 

X=2*PI •FREQ»RHO*CVEE 
X»X/FUDGE 
J0=MMBSJ0(X, IER) 

J1=MMBSJ1(X, IER) 

CALL  MMBSVN(X.0.0, l.YO.IER) 

CALL  MMBSVN(X, 0.999. 1, VI. IER) 

H0*CMPLX( J0,-V0( 1) ) 

H1*CMPLX( J1.-V1( 1) ) 

HR»H0/HJ. 

V»1 .0-X»Hfl 
V*X/V 

Y=Y»PI • 120*CVE£ 

V=-I»V 
VA=CABS( V ) 

PH= AT AN2 ( A  I  MAGI  V ) , RE AL ( V ) ) 

PH=PH» 180 . 00/PI 

WRITE (10,600) RHO . VA . PH 

RH0=RH0*INC 

IF(RHO.lE.MAXOIS)  GO  TO  100 

•FORMAT  STATEMENTS* 

F0RMAT(4X. 'NIGHT  (N)  OR  DAV  (0)7: ') 

FORMAT ( A2 ) 

FORMAT (9X. 'TOO  *  DAVTIME') 

F0RMAT(9X. 'TOO  =  NIGHT) 

F0RMAT(9X , ' E/H  FAR  FIELD  APPROXIMATION') 

FORMAT(5X, ' ENTER  ORIGINAL  DISTANCE  FROM  SOURCE:') 
FORMAT( F9 . 2) 

FORMAT(5X, 'ENTER  INCREMENT  FOR  THIS  DISTANCE:  ) 
F0RMATC5X ENTER  MAX  DISTANCE  TO  COMPUTE:') 

FORMATt 2X ENTER  FREQUENCY:') 

FORMAT( F7 . 2 ) 

F0RMAT(9X. 'HORIZONTAL  ELECTRIC  DIPOLE  FREQ  =  .F7.2) 

FORMAT ( IX. '  ' ) 

FORMAT ( 12X, 'DISTANCE' , 14X, 'MAGNITUDE' . 15X, 'PHASE' ) 
FORMAT ( 10X.E12.5, 10X.E12.5,  10X.E12.5) 


nnnonn  o  o  nooooooooo 


CREATED: 

LAST  UPDATE: 

BY: 

PURPOSE: 


21FEB8S 
07MAR8S 
A.  KUZEL 

THIS  PROGRAM  CALCULATES  THE  VALUE  OF  E/M 
OF  AN  ELF  WAVEFORM  OF  ONE  OF  A  SET  OF 
FREQUENCIES  AT  VARIOUS  DISTANCES. 

>=  200  KM,  FROM  THE  VEO  SOURCE.  THIS  IS 
THE  FAR  FIELD  APPROXIMATION. 


•DECLARE  VARIABLES* 


COMPLEX 

REAL 

real 

CHARACTER 


I .V, HO. HI. HR 

JO. J1,V0( 1) ,V1(1) , MMBS JO , MMBS J 1 , MAXDI S 
RHO.INC.PI , CVEE. FREQ. FUDGE. X.VA, PH 
TOO 


•INITIALIZE  CONSTANTS* 


PI i4* ATAN ( 1.0) 
I=CMPLX(0.0, 1.0) 


THE  FOLLOWING  SECTION  IS  the  INTERACTIVE  SECTION  OF  ThE  PROGRAM 
OF  The  PROGRAM.  THE  USER  CHOOSES  STARTING  DISTANCE.  INCREMENT.  AND 
ENOING  DISTANCE  TO  CALCULATE,  NIGHT  OR  DAV  VALUES.  AND  THE  FREQUENCY  . 
OF  INTEREST. 


WRITE(S.540) 

REAO(5,545) RHO 

WR I TE ( 5 . 547 ) 

READ ( 5 . S60 )  I NC 

FUOGE=3.0E*05 

WRITE(5,548) 

REAOt  5 . 545  JMAXOI S 
WRITE  I  5. 550) 

REAO(5.560)FREQ 
*R : 'E (  5 . 205  I 
READ! 5 , 206 ) TOD 

C  ’SELECT  PROPER  IONOSPHERIC  REFLECTION  HEIGHT  AND 

C/V  VALUES  BASED  ON  USER'S  INPUT  DATA* 

I F ( FREQ . EQ . 30 )  THEN 

I F ( TOO . EQ .  D' )  Then 
RFlhT  =  46 .  1 
CvEE= 1 . 34 
END  IF 

I F (TOO . EQ . ' N ' )  THEN 
RFLHT=72 . 0 
CvEE=  1  .  12 
END  IF 
ENO  IF 

I F / FREQ . EQ . 50 >  Then 
I F I '00  EQ .  0  ;  Then 

SF^-t-47  3 

C  vE £  =  l  . 30 

ENO  IF 

IF ( TOO . EQ . ' N  )  THEN 
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RFLHT =73.3 
CVEE*1. 11 
END  IF 
END  IF 

I F (FREQ . EO . 75 )  THEN 
I F (TOD . EQ . ' D ' )  THEN 
RFLHT  =  49 . 1 
CVEE* 1 . 27 
END  IF 

IF ( TOO . EQ . 'N ' )  THEN 
RFLHT*74 . 3 
CV£E*1. 10 
END  IF 
ENO  IF 

IF(FREQ.EQ. 100)  THEN 
IF(TOO . EQ . '0 ' )  THEN 
RFLHT*5Q . 1 
CVEE* 1 . 25 
END  IF 

I F (TOD . EO .  N' )  THEN 
RFLHT  =  75.0 
CVEE=1.09 
END  IF 
ENO  IF 

IF( FREQ. EO. 150)  THEN 
I  F  (TOD  .  EQ  .  '  0  '  )  THEN 
RFLHT  =51.4 
CVEE* 1.22 
ENO  IF 

I F ( TOD . EO . 'N' )  THEN 
RFLHT=76.0 
CVEE* 1.09 
ENO  IF 
ENO  IF 

I F ( FREQ . EQ . 200 )  THEN 
I F ( TOO . EQ . 'O' )  THEN 
RFLHT  =  52 . 4 
C vEE  = 1 . 20 
ENO  I* 

I  F  (  TOO  .  EQ  .  N  )  Then 
R  F  LHT  =  76 . 8 
CvEE* 1 . 08 
ENO  IF 
END  IF 

I F ( FREQ . EQ . 300 )  THEN 
IF  (  TOO  .  EQ  .  '  D  '  )  THEN 
RFLHT=53 . 7 
CvEE  =  1 .  18 
ENO  IF 

I F ( TOD . EQ . 'N ’ )  THEN 
RFLHT  =  77 . 8 
CvEE* 1.07 
ENO  IF 
ENO  IF 

IF(FREQ.EQ.400)  CvEE=1.  15 
IF (FREO  EO  800)  CvEE  =  1  1  1 
I F ( FREQ . EQ . 1600)  CvEE=1  07 

•write  headers  for  output  table* 


WRITE!  10. 220) 


TR  7521 


v  w  w  _  if,  ^  :  ~  ^  -  -  -  v  j  —  .  _  -  v- 


V  V  V  V 


MR  I TE ( 10 . S7S ) 

MRITEC 10,570)FR£Q 
I F  ( TOO . EQ .  ' 0 ' )  MRITEt 10.211) 

I F ( TOD . EQ . ' N ' )  MRITEt 10 . 212) 
MR  I  TE  t  10.575) 

MRITE( 10,580! 

MRITEt  10.575) 

C  •PERFORM  THE  CALCULATION* 

100  X»2*PI*FREQ*RH0*CVEE 

X»X/ FUDGE 
J0*MMBSJ0(X, !ER) 

J 1*MMBS J 1 ( X . IER) 

CALL  MMBSYNt  X , 0.0,1, VO. !ER) 
CALL  MMBSVNtX. 0.999, 1, VI. IER) 
H0*CMPLX( JO. -V0( 1)  ) 

H1=CMPLX( J1 . -Vl(  1)  ) 

HR=H0/M1 
V  =  HR 

VfV*PI*120*''VEE 
V=-I *V 
VA=CABS(V) 

PH=ATAN2(AIMAG(V  ) ,REAL( V)  ) 
PH=PH* 180. 00/PI 
MRITEt 1Q.600)RH0. VA.PH 
RHO=RHO* INC 

IF(RHO.LE.MAXOIS)  GO  TO  100 
C  ‘FORMAT  STATEMENTS* 


205 

FORMAT ( 4X , ' NIGHT 

tN)  OR  DAV  t  D) 7  ;  •  ) 

206 

FORMAT  t  A2 ) 

211 

FORMAT 1 9X ,  'TOO  * 

OAVTIME ' ) 

212 

F0RMAT(9X. 'TOO  = 

NIGHT- ) 

220 

FQRMAT(9X, • E/H  FAR  FIELD  APPROXIMATION’) 

540 

FORMAT ( 5X . • ENTER 

DISTANCE  from  SOURCE:) 

545 

FORMAT! F9 . 2 ) 

547 

FORMA  T ( 5x .  • ENTER 

INCREMENT  FOR  DISTANCE:  ) 

546 

F0RMAT(5X. 'ENTER 

MAXIMUM  distance  TO  SE  COMPUTED: 

550 

FORMATt  2X ,  ENTER 

FREQUENCY : ' ) 

560 

FORMAT (F7. 2) 

570 

F0RMAT(9X, ’VERTICAL  ELECTRIC  DIPOLE  FREQ 

=  ’  .  F7 

575 

FORMATt IX.' 

) 

580 

FORMATt 12X , 'DISTANCE ' . 14X , ' MAGNITUDE ’ , 14X . 

PHASE 

600 

FORMATt 10X.E12.5 

.  I0X.E12.5, 10X.E12.5) 

ENO 


ooooooo  r>  o  o  o  n  o  o  ooononn 


TR  7521 


CREATED: 

LAST  UPOATE: 
8V: 

PURPOSE: 


27FEB85 
27FEB85 
A.  KUZEL 

FORTRAN  SUBROUTINE  TO  CALCULATE  THE 
HYPERBOLIC  COTANGENT  OF  THE  PASSED  PARAMETER. 


SUBROUTINE  COTH(ZZ) 

REAL  ZZ.XX 

XX=EXP(ZZ)*EXP<-ZZ) 

XX*XX/(EXP(ZZ)-EXP(-ZZ) ) 

ZZ*XX 

RETURN 

END 


CREATED: 

LAST  UPOATE: 
BY: 

PURPOSE: 


27FEB85 
27FEB8S 
A.  KUZEL 

FORTRAN  SUBROUTINE  TO  CALCULATE  THE 
HYPERBOLIC  COSECANT  OF  THE  PASSED  PARAMETER. 


SUBROUTINE  CSCH(ZZ) 

REAL  ZZ.XX 

XX  =  2 

XX=XX/(EXP(ZZ)-EXP(-ZZ) ) 

ZZ  =  XX 

RETURN 

END 


CREATED: 

LAST  UPDATE: 
BY: 

PURPOSE: 


27FEB85 
27FEB85 
A.  KUZEL 

FORTRAN  SUBROUTINE  TO  CALCULATE  THE 
HYPERBOLIC  TANGENT  OF  THE  PASSED  PARAMETER. 


SUBROUTINE  TANH(ZZ) 

real  ZZ.XX 

XX=EXP(ZZ)-EXP(-ZZ) 
XX=XX/(EXP(ZZ)*EXP( -ZZ) ) 

zz=xx 

RETURN 

END 
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